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Abstract. The quark-hadron duality is studied in a systematic way for polarized and unpolarized structure 
functions, by taking into account all the available data in the resonance region. In both cases, a precise 
perturbative QCD based analysis to the integrals of the structure functions in the resonance region has 
been done: non perturbative contributions have been disentangled and the higher twist contributions have 
been evaluated. A different behavior for the unpolarized and polarized structure functions at low Q 2 has 
been found. 
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1 Introduction 

Understanding the structure and interaction of hadron in 
terms of the quark and gluon degrees of freedom of QCD 
is one of the unsolved problems of the Standard Model 
of nuclear and particle physics. At present it's not possi- 
ble to describe the physics of hadrons directly from QCD, 
however it is known that it should just be a matter of 
convenience the choice of describing a process in terms of 
quark-gluon or hadronic degrees of freedom. This concept 
is called quark-hadron duality. At high energies, where the 
interactions between quarks and gluons become weak and 
quark can be considered asymptotically free, an efficient 
description of phenomena is possible in terms of quarks. 
At low energies, where the effects of confinement become 
large, it is more efficient to work in terms of collective 
degrees of freedom, the physical mesons and baryons. In 
these terms, it's clear that the duality between the quark 
and hadron descriptions reflects the relationship between 
confinement and asymptotic freedom, and is intimately re- 
lated to the nature of the transition from non-perturbative 
(low energy) to perturbative QCD (high energy). 

The concept of duality was introduced for the first 
time by Bloom and Gilman [T] in deep inelastic scattering 
(DIS). They noticed an equivalence between the smooth 
x dependence of the inclusive structure function at large 
Q 2 and the average over W 2 of the nucleon resonances. 
Furthermore, this equivalence appeared to hold for each 
resonance, over restricted regions in W. Based on this ob- 
servations, one can refer to global duality if the average, 
defined e.g. as the integral of the structure functions, is 
taken over the whole resonance region 1 < W 2 < 4 GeV 2 . 
If, however, the averaging is performed over smaller W 2 
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ranges, extending e.g. over single resonances, one can refer 
to local duality. 

Although the duality between quark and hadron de- 
scriptions is, in principle, formally exact, how this reveals 
itself specifically in different physical processes and under 
different kinematical conditions is the key to understand 
the consequences of QCD for hadronic structure. The phe- 
nomenon of duality is quite general in nature and can be 
studied in a variety of processes, such as DIS, e + e _ an- 
nihilation into hadrons, and hadron-hadron collisions, or 
semi-leptonic decays of heavy quarks. With the advent of 
both more detailed studies of soft scales and confinement 
[2], and higher precision measurements covering a wide 
range of reactions, it is now becoming possible to investi- 
gate the role of duality in QCD as a subject per se. 



2 Kinematical variables 

Besides the scaling variable x, other variables have been 
used in the literature to study duality and a number of 
parametrizations based on these variables have been pro- 
posed that reproduce in an effective way some of the cor- 
rections to the perturbative QCD calculations. The most 
extensively used variables are: x' — where to' = 

l/x+M 2 /Q 2 , originally introduced by Bloom and Gilman 
in order to obtain a better agreement between DIS and 
the resonance region; £ = 2x/(l + (1 + Ax 2 M 2 /Q 2 ) 1 / 2 ) 
originally introduced to take into account the target mass 
effects; x w = Q 2 + B/(Q 2 + W 2 - M 2 + A), A and B be- 
ing fitted parameters, used in Refs. 01121 • These additional 
variables include a Q 2 dependence that phenomenologi- 
cally absorbs some of the scaling violations that are impor- 
tant at low Q 2 . In Fig^Jtheir behavior vs. x is compared 
for different values of Q 2 . From the figure one can see that 
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Fig. 1. Ratio between the three variables x', £ and Xw defined 
in the text and the Bjorken variable i as a function of x. 



the numerator is evaluated using the experimental data in 
the resonance region, while the one at the denominator is 
calculated from parametrizations that reproduce the DIS 
behavior of the data at large Q 2 . The ratios have been 
calculated in unpolarized and polarized cases. It has been 
found ^21 that quark-hadron duality has not been ful- 
filled by using solely the parton distribution functions up 
to NLO in both the unpolarized and polarized structure 
functions F 2 and g\ . However it was possible to see a differ- 
ent behavior between i? U npoi and i? po i- In the unpolarized 
case the ratio is increasing with Q 2 , but for the polarized 
case the situation is different: while at low Q 2 the ratio 
is significantly below unity and shows a strong increase 
with Q 2 , at higher Q 2 the ratio derived from HERMES is 
above unity and it appears to be weakly dependent of Q 2 
within error bars. The situation is different with the use 
of the phenomenological fits to DIS data 14,15,5 . Since 
these phenomenological parametrizations are obtained by 
fitting deep-inelastic data even in the low Q 2 region, they 
can implicitly include non-perturbative effects and this 
may explain the "observation of duality". It becomes re- 
ally important to understand the contribution of these 
non-perturbative effects. 



by calculating F2 in £ and x' , one effectively "rescales" the 
structure function to lower values of x, in a Q 2 dependent 
way, namely the rescaling is larger at lower Q 2 . In the 
analysis reported in the following, the x variable has been 
used and all the corrections have been applied one by one. 



3 Analysis of data 

A quantitative analysis of the Q 2 dependence of quark- 
hadron duality in both polarized and unpolarized ep scat- 
tering is presented. All current data in the resonance re- 
gion, 1 < W 2 < 4 GeV 2 , have been taken into account. 
For the unpolarized case it has been used the data ob- 
tained at Jefferson Lab in the range 0.3 < Q 2 < 5 GeV 2 
UJ, and the data from SLAC ([7] and references therein) 
for Q 2 > 4 GeV 2 . For the polarized case there are only 
few experimental data in the resonance region. One set is 
part of the E143 data [5], and it corresponds to Q 2 = 0.5 
and 1.2 GeV 2 . Another set is the one from HERMES |3 
□HI in the range 1.2 < Q 2 < 12 GeV 2 . 

In the polarized case the Q 2 dependence originates 
from the structure function F\ and from the ratio R, for 
which the SLAC global analysis parametrization has 
been used. For the asymmetry A±, it was used a power 
law fit to the world DIS data at x >0.3, A\ — x°- 7 , as 
already shown in Ref. This parametrization of A\ is 
constrained to 1 at x—1 and it does not depend on Q 2 , as 
indicated by experimental data in this range |12| . 

The full procedure of the analysis is described in ^3] ■ 
The quark-hadron duality in DIS is studied by consider- 
ing the ratio of the integrals of the structure functions 
integrated in a defined ir-range, corresponding to the W 
range of the resonance region. The structure function in 



4 Size of non-perturbative contributions 

In order to understand the nature of the remaining Q 2 
dependence that cannot be described by NLO pQCD evo- 
lution, the effect of target mass corrections and large x 
resummation have been studied. As mentioned early, the 
analysis was performed by using x as an integration vari- 
able, which avoids the ambiguities associated to the us- 
age of other ad hoc kinematical variables. Standard in- 
put parametrizations with initial scale Q 2 = 1 GeV 2 have 
been used. Once both effects have been subtracted from 
the data, and assuming the validity of the twist expansion, 
one can interpret any remaining discrepancy of the ratio 
from unity in terms of higher twist. 

4.1 Target Mass Corrections (TMC) 

TMC are necessary to take into account the finite mass 
of the initial nucleon. They are corrections to the leading 
twist (LT) part of the unpolarized structure function F2. 
For Q 2 larger than « 1 GeV 2 , TMC are taken into account 
through the following expansion |16| : 

^ MC (-^ 2 ) = i $^,Q 2 ) + 6^^p (1) 

where F%° is the structure function in the absence of TMC. 
Following the original suggestion of fT7| , only terms up to 
order M 2 /Q 2 are kept in the expansion, so as to minimize 
ambiguities in the behavior of F 2 at x » 1. Although this 
procedure disregards parton off-shell effects that might be 
important in the resonance region flSHI)] . it's important 
to emphasize its power expansion character, and setting 
as a limiting condition for its validity, that the inequality 
x 2 M 2 /Q 2 < 1 be verified [20], Q 2 ^ 1.5 GeV 2 . 



Alessandra Fantoni, Nicola Bianchi, Simonetta Liuti: Quark-Hadron Duality in Structure Functions 



3 



4.2 Large x Resummation (LxR) 

LxR effects arise formally from terms containing powers 
of ln(l — z), z being the longitudinal variable in the evolu- 
tion equations, that are present in the Wilson coefficient 
functions C'(z). The latter relate the parton distributions 
to e.g. the structure function F2, according to: 



H 



F 2 NS (x, Q 2 ) = g jf dz C NS (z) q NS (x/z, Q 2 ), 



(2) 



where it has been considered only the non-singlet (NS) 
contribution to F2 since only valence quarks distributions 
are relevant in the present kinematics. The logarithmic 
terms in Cns(z) become very large at large x, and they 
need to be resummed to all orders in cvg. This can be ac- 
complished by noticing that the correct kinematical vari- 
able that determines the phase space for the radiation of 
gluons at large x, is W 2 — Q 2 (l — z)/z, instead of Q 2 (2U 
I22j . As a result, the argument of the strong coupling con- 
stant becomes z-dependent: as(Q 2 ) — > as(Q 2 (l — z)jz) 
( |23| and references therein). In this procedure, however, 
an ambiguity is introduced, related to the need of con- 
tinuing the value of as for low values of its argument, 
i.e. for z very close to 1 [21] ■ The size of this ambiguity 
could be of the same order of the HT corrections. Nev- 
ertheless, the present evaluation is largely free from this 
problem because of the particular kinematical conditions 
in the resonance region. In this analysis, in fact, the struc- 
ture functions have been studied at fixed W 2 , in between 
1 < W 2 < 4 GeV 2 . Consequently Q 2 increases with x. 
This softens the ambiguity in as, and renders this proce- 
dure reliable for the extraction of HT terms. 



4.3 Disentangle of non-perturbative contributions 

All the effects described in the present section are sum- 
marized in Fig|2 where the ratio between the resonance 
region and the 'DIS' one is reported for the unpolarized 
and for the polarized case: the numerator is obtained from 
the experimental data, while the denominator includes the 
different components of the present analysis, one by one. 

For unpolarized scattering it has been found that TMC 
and LxR diminish considerably the space left for HT con- 
tributions. The contribution of TMC is large at the largest 
values of Q 2 because these correspond also to large x val- 
ues. Moreover, the effect of TMC is larger than the one of 
LxR. The lowest data point at Q 2 rs 0.4 GeV 2 has been ex- 
cluded from the analysis because of the high uncertainty in 
both the pQCD calculation and the subtraction of TMC. 

Similarly, in polarized scattering the inclusion of TMC 
and LxR decreases the ratio R^Ji- However, in this case 
these effects are included almost completely within the 
error bars. Clearly, duality is strongly violated at Q 2 < 1.7 
GeV 2 . 

The difference between unpolarized and polarized scat- 
tering at low Q 2 can be attributed e.g. to unmeasured, so 
far, Q 2 dependent effects, both in the asymmetry, A±, and 
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Fig. 2. Ratio between the integrals of the measured structure 
functions and the calculated ones plotted as a function of Q 2 . 
The calculation includes one by one the effects of NLO pQCD 
(squares), TMC (open circles) and LxR (triangles), The top 
panel refers to the unpolarized case, while the bottom panel to 
the polarized one. 



in 52 • Furthermore, a full treatment of the Q 2 dependence 
would require both a more accurate knowledge of the ratio 
R in the resonance region, and a simultaneous evaluation 
of 52 ■ The present mismatch between the unpolarized and 
polarized low Q 2 behavior might indicate that factoriza- 
tion is broken differently for the two processes, and that 
the universality of quark descriptions no longer holds. 
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5 Size of Higher Twist (HT) corrections 

The discrepancy from unity of the ratios already presented 
is interpreted in terms of HTs. In Figs. 13141 the question 
of the size of the HT corrections is addressed explicitely. 
For F 2 , they are defined as: 



H(x,Q 2 



Q 



(F™(x,Q 2 )-F™) 



Cht(x) 



H(x,Q 2 



7 pQCD 



(x/Q 2 



Q 
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2*2 



(x,Q 2 ) 



F 2 LT 



(4) 



A similar expression is assumed for g±. Cht is the so- 
called factorized form obtained by assuming that the Q 2 
dependences of the LT and of the HT parts are similar 
and therefore they cancel out in the ratio. Although the 
anomalous dimensions of the HT part could in principle 
be different, such a discrepancy has not been found so 
far in accurate analyses of DIS data. The HT coefficient, 
Cht has been evaluated for the three cases listed also in 
FigUl namely with respect to the NLO pQCD calculation, 
to NLO+TMC and to NLO+TMC+LxR. The values of 
1 + Crt/Q 2 are plotted in FigEI (upper panel) as a func- 
tion of the average value of x for each spectrum. One can 
see that the NLO+TMC+LxR analysis yields very small 
values for Cht in the whole range of x. Furthermore, the 
extracted values are consistent with the ones obtained in 
Ref. [2131 using a different method, however the present ex- 
traction method gives more accurate results. Because of 
the increased precision of our analysis, we are able to dis- 
entangle the different effects from both TMC and LxR. 

In the polarized case (lower panel) the HTs are small 
within the given precision, for Q 2 > 1.7 GeV 2 , but they 
appear to drop dramatically below zero for lower Q 2 val- 
ues. The inclusion of TMC and LxR renders these terms 
consistent with zero at the larger Q 2 values, but it does not 
modify substantially their behavior at lower Q 2 . It should 
also be noticed that, by parametrizing the structure func- 
tions as in Eq.@, it is assuming that all of the non- 
perturbative (np) contributions are included in 0(1/ Q 2 ) 
twist-4 terms. These are in fact the largest type of devia- 
tions from a pQCD behavior, to be expected at Q 2 values 
of the order of few GeV. Only from accurate analyses us- 
ing a larger number of more precise data, would one be 
able to distinguish among different np behaviors. From a 
comparison with results of ratio including phcnomcnologi- 
cal parametrizations |13| that includes some of these extra 
np behaviors it's possible to see, however, that their effect 
seems not be large. 

In Fig0]the results of the present analysis in the unpo- 
larized case are compared to other current extractions of 
the same quantity. These are: i) the extractions from DIS 
data, performed with the cut W 2 > 10 GeV 2 [55ll26ll27| : 
ii) the recent DIS evaluation by S. Alekhin [SS] using a cut 
on W 2 > 4 GeV 2 , and including both TMC and NNLO; 
in) the results obtained within a fixed W 2 framework |2U| . 
including both TMC and LxR. The results obtained in the 
deep inelastic region [21] also including both TMC and 
LxR yield small HT coefficients, consistent with the ones 
found in Ref. 20 . However, while most of the suppression 



of the HT in the resonance region is attributed to TMC, in 
[SHj the contribution of TMC is small and the suppression 
is dominated by LxR. In other words, the Q 2 behavior in 
the DIS and resonance regions seems to be dominated by 
different effects. 



(3) 6 Conclusions 



A precise and detailed analysis of all publish data in reso- 
nance region has been presented, with the aim of studying 
the quark-hadron duality in unpolarized and polarized ep 
scattering. A pQCD NLO analysis including target mass 
corrections and large x resummation effects was extended 
to the integrals of both unpolarized and polarized struc- 
ture functions in the resonance region. Both effects have 
been quantified and disentangled for the first time. In the 
present analysis ^3], duality is satisfied if the pQCD cal- 
culations agree with the data, modulo higher twist con- 
tributions consistent with the twist expansion. A different 
behavior for unpolarized and polarized structure functions 
has been found, and duality seem strongly violated in the 
latter case for Q 2 <1.7 GeV 2 The discrepancy of the ratio 
from unity has been interpreted in terms of HTs. While 
the size of the HT contributions is comparable in both 
polarized and unpolarized scattering at larger x and Q 2 
values, at low x and Q 2 large negative non-perturbative 
contributions have been found only in the polarized case. 
The present detailed extraction of both the Q 2 depen- 
dence and the HTs in the resonance region establishes a 
background for understanding the transition between par- 
tonic and hadronic degrees of freedom. In particular, it 
seems to be detecting a region where the twist expansion 
breaks down, and at the same time, the data seem to be 
still far from the Q 2 — ► limit, where theoretical predic- 
tions can be made [HOj- More studies addressing this re- 
gion will be pursued in the future, some of which are also 
mentioned in |20II31| . A breakdown of the twist expansion 
can be interpreted in terms of the dominance of multi- 
parton configurations over single parton contributions in 
the scattering process. In order to confirm this picture it 
will be necessary to both extend the studies of the twist 
expansion, including the possible Q 2 dependence of the 
HT coefficients [321133] and terms of order 0(1/Q 4 ), and 
to perform duality studies in semi-inclusive experiments. 
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